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Simulations in the e-Regime of Chiral Perturbation Theory * 
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We discuss the potential of Ginsparg- Wilson fermion simulations in the e-regime of chiral perturbation theory, 
regarding the determination of the leading low energy constants of the effective chiral Lagrangian. It turns out to 
be very hard to measure observables in the topologically trivial sector. There a huge statistics would be required, 
due to the frequent occurrence of very small eigenvalues. Moreover, contact with chiral perturbation theory is 
only established if the physical volume of the system is sufficiently large. 
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QCD at low energy can be described by an ef- 
fective Lagrangian C g for the quasi- Goldstone 
bosons of chiral symmetry breaking. £ c g may 
then be evaluated by chiral perturbation theory 
(xPT). The case of a finite volume and very light 
quarks — so that the pion Compton wave length 
clearly exceeds the box length — is called the 
e-regime. There the zero-mode is treated by col- 
lective variables, while the excitations (which do 
fit into the box) can be evaluated perturbatively 
in the e-expansion [1]. In this framework, ob- 
servables strongly depend on the topology of the 
gauge background [2] , hence they should be mea- 
sured in fixed topological sectors. 

£ e ff involves coupling constants which appear 
as free parameters in x?T. In principle they 
can be determined by comparison with numerical 
data from lattice QCD in the e-regime. The same 
values are also relevant for the physical situation 
in a large volume (the standard \PT), hence their 
knowledge is very important in view of experi- 
ments and further numerical studies, especially 
for the chiral extrapolation. Their numerical eval- 
uation, however, may be easier in the e-regime, 
due to the usage of a small volume. 

Such simulations in the e-regime are now 
conceivable by using Ginsparg- Wilson fermions, 
which avoid additive mass renormalization and 
define the topological charge through the Index 
Theorem [3]. Here we discuss the question how 
one could extract results for the low energy con- 
stants which occur in the leading order of C c g, the 



pion decay constant F v and the scalar condensate 
£. 

We performed quenched QCD simulations with 
the Wilson gauge action and the overlap Dirac 
operator [4] 



D m =(l-£)D<$+m, 

£>$ - M(l + A/VA^A) , A = D w - fi 



(1) 



*Poster presented by K.-I. Nagai at Lattice 2003 
Preprint DESY 03-126, HU-EP-03/56, SFB/CPP-03-31 



where m is the quark mass and Dw the Wil- 
son operator. At (3 — 6 and j3 = 5.85 we set 
jj, = 1.4 resp. /j = 1.6. We approximated the in- 
verse square root by Chebyshev polynomials to 
an accuracy of 10~ 12 . To identify the low ly- 
ing eigenvalues (EVs) of D$ ^Dov we applied 
alternatively the Ritz functional method and the 
Arnoldi algorithm. The latter allows for the de- 
termination of many EVs (0(100)). The very 
lowest EVs determine the index v. there is usu- 
ally a gap by several orders of magnitude between 
the numerical values of the zero EVs and the rest 
of the spectrum. The chirality of the zero modes 
reveals the sign of v. Two typical examples are 
given in Table 1. 

On a 10 3 x 24 lattice at j3 = 6 we also com- 
pared the index with the topological charge de- 
termined from (standard) cooling. In this respect, 
we investigated 51 configurations: for 41 of them 
the index of D ov and the charge determined from 
cooling coincided, and in the remaining 10 cases 
it differed by 1. 

Next we mention that Random Matrix Theory 
(RMT) has been applied to QCD and it yields 
predictions for the probability distributions of the 
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EVs oiDZ n D^ 


chirality 


example for charge v = 


6.23013e- 3 


0.3608201 


6.23014e- 3 


-0.3608198 


9.90174e- 3 


-0.7622260 


9.90177e- 3 


0.7622260 


2.69086e- 2 


-0.7528350 


example for char 


ge v = — 2 


6.03761e- 16 


-1.0000000 


9.43764e- 11 


-0.9999999 


3.68203e - 3 


-0.9240370 


3.73623e- 3 


0.9205659 


7.04176e - 3 


-0.7016518 



Table 1 

Typical examples for the lowest eigenvalues of 
£>ov f £>ov on a 10 3 x 24 lattice at (3 = 6. 

low lying eigenvalues [5]. For the first non-zero 
EV these predictions are confirmed if the physi- 
cal volume is sufficiently large, V>(1.2 fm) 4 [6]. 
In particular, we know from these probability dis- 
tributions that in the topologically trivial sector 
there is a significant density of very small (non- 
zero) EVs. As an example, we found a config- 
uration on the 10 3 x 24 lattice at j3 = 6 where 
the lowest EVs of D$ read: 7.32e - 5, 

7.33e — 5, 1.04e — 4 etc. The corresponding chi- 
ralities -0.127, 0.127, 0.826 show clearly that 
this configuration has index 0. This is also con- 
firmed from cooling, but the form of the cooled 
configuration suggests that charge arises from 
the cancellation of a pair of topological objects 
with charges ±1, see Fig. 1. The index never in- 
volved such a cancellation. However, in the case 
of small non-zero EVs there seems to be a trend 
towards a cancellation picture from cooling, as in 
the example discussed here. 

We now address the meson correlators as phys- 
ical observables, from which one can try to ex- 
tract F n and E. The corresponding formulae for 
quenched have been worked out in Ref. [7]. 
For a first comparison to numerical data we refer 
to Refs. [8]. The axial current correlation func- 
tion on an L 3 x T lattice, 



(A M (t)A M (0))„ = 
2hi(r) = r 2 -r 
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^ + 2m £ M CO T-/u(t) 
1/6, r = t/T 



(2) 



Figure 1. A 2-d cut of the configuration discussed 
in the text after 200 cooling sweeps. The vertical 
axis shows the naive topological charge density, 
which extends over both signs, suggesting some- 
thing like an instanton anti-instanton pair. 

^-=z[l v {z)K v {z)+I v+1 {z)K v _ 1 {z) 

is most suitable for our purpose, because it only 
involves F v and £ as free parameters. The deter- 
mination of F n is relatively easy in the minimum 
at t = 1/2. E is related to the curvature and far 
more difficult to extract, since the sensitivity of 
this curvature to a change of E is hardly visible 
over a wide range 2 . However, we saw also here 
that the physical volume must exceed the lower 
limit given before in the context of RMT, other- 
wise one does not obtain the curvature required 
by any acceptable value of E. Once we are in the 
right regime, v = may look like the simplest 
case from eq. (2), but from the numerical point 
of view this is a nightmare. The problems in this 
sector are related to the danger of very small EVs. 
Ref. [9] pointed out before that they may be a se- 
rious problem, and our observations confirm this 
in a striking way. As an example, we show his- 
tories for the measurement of {A (t) A (0))o and 
(A (t) A)(0))i in Fig. 2. In the topologically neu- 
tral sector there are strong spikes at those config- 
urations which have very small EVs. Their height 
is maximal at very small quark mass, as Fig. 3 il- 
lustrates. 

Obviously the occurrence of such spikes makes 
it very hard to measure expectation values - 
they stabilize only with a huge statistics. To get 



2 This situation is somehow complementary to the RMT 
comparison of the EVs, where E is easier to determine [6]. 
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Figure 2. Histories of the axial correlation func- 
tion Ao(t)Ao(0) at t = 8 and at t = 16, with 
m = 21.3 MeV. For some v = configura- 
tions (on the left) we recognize pronounced spikes. 
At \v\ = 1 (on the right) the history is much 
smoother. 



8e-4 



6e-4 



4e-4 



2e-4 



- conf. 20 
« conf. 12 



5 10 15 20 

m [MeV] 

Figure 3. The height of the two spikes of 
A (16)^4o(0) in the v = history of Fig. 2, at 
different quark masses m. 

an estimate for this effect, we consider the con- 
tribution of the smallest (non-zero) EV alone to 
the scalar condensate. This contribution reads 

2m 



?%L = h r dzP y {z)- 
v Jo 1 



where we insert the probabilities P v {z) provided 
by RMT. We now generated fake values for the 
lowest EV with probability P u (z) and computed 
EwjL from them. The result can only be trusted 
when the standard deviation becomes practically 
constant. Fig. 4 shows the evolution of this stan- 
dard deviation as the statistics is enhanced. We 
see that the sector v = would require a tremen- 
dous statistics of O(10 4 ) configurations, but the 
situation is clearly better for v =^ 0. 
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Figure 4. The standard deviation in the (fake) 
measurement ofTr")„, the contribution of the low- 
est non-zero EV to the scalar condensate. 

We conclude that there are several conditions 
for running conclusive simulations with Ginsparg- 
Wilson fermions in the e-regime. The volume 
should obey V>(1.2 fm) 4 , and measurements 
have to be performed in a sector of non-trivial 
topology. The quark mass should be small for 
conceptual reasons, but taking it too small causes 
technical problems in the measurements. 
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